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We consider the polarization arising from scattering in an en- 
velope illuminated by a central anisotropic source. This work ex- 
tends the theory introduced in a previous paper (Al-Malki et al. 
1999) in which scattering polarization from a spherically sym- 
metric envelope illuminated by an anisotropic point source was 
considered. Here we generalize to account for the more realis- 
tic expectation of a non-spherical envelope shape. Spherical har- 
monics are used to describe both the light source anisotropy and 
the envelope density distribution functions of the scattering par- 
ticles. This framework demonstrates how the net resultant polar- 
ization arises from a superposition of three basic "shape" func- 
tions: the distribution of source illumination, the distribution of 
envelope scatterers, and the phase function for dipole scattering. 
Specific expressions for the Stokes parameters and scattered flux 
are derived for the case of an ellipsoidal light source inside an 
ellipsoidal envelope, with principal axes that are generally not 
aligned. Two illustrative examples are considered: (a) axisym- 
metric mass loss from a rapidly rotating star, such as may apply 
to some Luminous Blue Variables, and (b) a Roche-lobe filling 
star in a binary system with a circumstellar envelope. As a gen- 
eral conclusion, the combination of source anisotropy with dis- 
torted scattering envelopes leads to more complex polarimetric 
behavior such that the source characteristics should be carefully 
considered when interpreting polarimetric data. 

Polarization - Stars: circumstellar - Stars: binaries - Stars: 
rotation - Stars: winds 



1. Introduction 

For many stars linear polarization is produced mainly from scat- 
tering of starlight by circumstellar matter (Rruszewski et al. 
1968; Serkowski 1970; Dyck et al. 1971; Shawl 1975). This po- 
larization can be used as a diagnostic of the geometry of the cir- 
cumstellar envelope and of the light source(s) (e.g., Shakhovskoi 
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1965; Serkowski 1970; Brown & McLean 1977; Brown et al. 
1978; Rudy & Kemp 1978; Simmons 1982, 1983; Friend & 
Cassinelli 1986; Clarke & McGale 1986, 1987). In many models 
for circumstellar scattering polarization, the illuminating sources 
are treated as isotropic point sources (e.g., Brown & McLean 
1977; Brown et al. 1978; Rudy & Kemp 1978; Shawl 1975; 
Simmons 1982). The effect of the finite size of the star as a 
light source has been studied (Cassinelli, Nordsieck, & Murison 
1987), including the effects of limb darkening (Brown et al. 
1989) and stellar occultation (Brown & Fox 1989; Fox & Brown 
1991; Fox 1991). Variable polarization is also revealing about a 
system, as evidenced in a recent study by Elias, Koch, & Pfeiffer 
(2008). With recent emphasis on clumped wind flows of hot 
stars (e.g., Hamann, Feldmeier, & Oskinova 2008), models have 
been developed to interpret variable polarization from hot star 
winds (Richardson, Brown, & Simmons 1996; Brown, Ignace, 
& Cassinelli 2000; Li et al. 2000; Davies et al. 2007). 

Relatively little has been done to explore the anisotropy of 
the source illumination and its consequences for interpreting 
polarimetric observations. There are some exceptions, such as 
a study of gravity-darkening effects for the polarization from 
Be stars by Bjorkman & Bjorkman (1994), the influence of star 
spots for polarizations from pre-main sequence stars by Vink 
et al. (2005), and more recently calculations of source anisotropy 
for interpreting polarimetric behavior observed in the post-red 
supergiant star HD 179821 by Patel et al. (2008). In a previ- 
ous paper (Al-Malki et al. 1999; hereafter Paper I), we modeled 
the polarization arising from Thomson and Rayleigh scattering 
explicitly for an arbitrary anisotropic (point) light source. As a 
proof of concept, the circumstellar envelope was taken as spher- 
ically symmetric to explore the polarization signals that arise 
solely from the properties of the source. An upper limit to the 
polarization of about 20% was derived in the limiting case of a 
disk-like star when viewed edge-on. 

This paper generalizes the results of Paper I by allowing both 
the point light source and the particle density distribution to be 
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functions of arbitrary angular form. As in Paper I, we restrict 
ourselves to optically thin envelopes and single Thomson or 
Rayleigh scattering mechanisms. Although the effects of a finite- 
sized illuminator will be important for obtaining quantitatively a 
more accurate polarization amplitude (e.g., Ignace, Henson, & 
Carson 2008), the point source approximation is certainly ade- 
quate for exploring polarization trends of generalized source and 
envelope geometries. In Sect. 2, we develop the theory required 
to evaluate the polarization, based on Paper I and the discus- 
sion of Simmons (1982). Applications to specific cases of astro- 
physical interest are explored in Sect. 3, with emphasis given to 
a treatment of the star and the circumstellar envelope as ellip- 
soidal. A discussion of these results and the need for ongoing 
modeling is given in Sect. 4. 



2. General expression for scattered flux and 
polarization 

As in Paper I, we neglect the effects of finite star depolarization 
and stellar occultation. Consequently scattering particles in the 
circumstellar envelope "see" a point star. However, we model the 
anisotropy of the illumination with a "shape" function to repre- 
sent directional flux. 

Under these approximations, our goal is to derive the net 
polarization from dipole scattering for an unresolved system 
that allows for an arbitrary circumstellar envelope geometry and 
source geometry. Our approach is to represent the different ge- 
ometrical factors in terms of expansions in spherical harmonics. 
As a result, it is important to be clear about the different coordi- 
nate systems being employed. 

Following the notation of Paper I, a description of the prob- 
lem requires three coordinate systems, each centered on the star: 
one for the observer, one for the envelope, and one for the illumi- 
nation pattern of the star itself. In general, none of these systems 
are collinear. The systems are defined as follows: 

i) For the observer reference frame we define Cartesian coordi- 
nates (x, y, z) centered at the star, with spherical coordinates 
(r, 8, (p), the line-of-sight being the Oz-axis. Thus z is a unit 
vector in the direction of the observer from the star and x,y 
are observer coordinates in the plane of the sky. For a scat- 
tering point in direction f, the scattering angle is given by 
cos6> = z ■ r, as in Paper I. The angle <p - tan^y/x is the 
observer's polarization angle (orientation) for any scattering 
point. 

ii) The star's frame is (X, Y, Z) with spherical coordinates 
(r, §, ip), where OZ is a convenient stellar axis (such as a ro- 
tation axis) that lies in the x—z plane of the observer's frame. 
The star system is rotated relative to the observer coordinates 
(x, y, z) through standard Euler angles (a, f}, y). 

iii) The envelope frame (X',Y',Z') with spherical coordinates 
(r,&',ip'), centered on the star, where OZ' is a convenient 
axis for the envelope. 

In general the scatterer density is n(r, <p'), and the flux of 
radiation from the star is F(r, §, tp). The stellar illumination is 
taken to be unpolarized. In fact, stellar atmospheres may show 
some low level of intrinsic polarization (Collins 1970); however, 
the polarization from the circumstellar scattering will be dom- 
inated by the Stokes-I component from the source. Following 
Eq. (1) from Paper I (see also Simmons 1982), the scattered flux 
and Stokes parameters (F sc , Q, U) of the scattered radiation at 
the Earth (distance D) are given in the form: 
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tp) r 2 x 



0'i + h) 

(h - h) exp(-2i'0) 



dr sin 8 d0 dip, 



(1) 



with the same definitions as in Paper I, with Z* — Q — iU, for 
i = V 11 !, k = 2n/A the wave number, and i\ and ii the scat- 
tering functions as defined by van de Hulst (1957). The circular 
polarization V is assumed zero. For Thomson (free electrons) or 
Rayleigh scattering, we have 



3 k 2 

h ±h = — cr(l ±cos 2 6»), 



(2) 



where the value of the cross section factor cr is chosen according 
to whether Thomson scattering (cr independent of k) or Rayleigh 
scattering (cr oc k 4 ) is considered. 

Provided that F varies smoothly, it may be expressed in 
terms of spherical harmonics in the observer frame (see Paper I): 

F(r, = ^ F lm {r) J] K m (a,/3, y) Y ln (8, <p), (3) 



1=0 m=-/ 



n=-/ 



where 

J F(r,O,<t>)Yj m (e,<p)d(cos0)d<f>, 



(4) 



and R l nm represents elements of a rotational matrix (see Paper I). 

Similarly, the density distribution of scatterers can also be 
expanded in terms of spherical harmonics, which becomes 



n'=l' 



n(r,e,<fi) = Yj n l'm'(r)Yi w (6,(/>), 
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where 



l'=0 m!=-V 



XI pin 
J n(r,6',ip')Y; m ,(0',(f>')d(cosO')d(f>'. (6) 



The Yi n (8, ip) factors in Eqs. (f3]l and (O are as defined in Paper I 
(c.f., Jackson 1975). For the product of two spherical harmonics, 
we can express the multipoles of n(r, 8, (p) and F(r, 8, <p) as: 



n(r,8,</>)F(r,8,4>) = n,, m ,(r) F, m (r) X 

IV mm' 

Y,, m ,(8,cp)Y Im (6,<p), 

or as 

n(r, 8, cp) F(r, M> = £ < m 2 C^ m , Y LM {8, 0). 

Imn I'm' 



(7) 
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In this latter expression, the set of C^ m , values are Clebsh- 
Gordon coefficients, arising from the products of two spherical 
harmonics. Only terms satisfying the following two conditions 
contribute to the sum in Eq. ((SJ: 



n + m' = M, 



(9) 
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Fig. 1. Illustration of the ellipsoidal envelope used in our models. The scattering envelope is shown in cross-section as shaded, 
with an equatorial width ranging from an inner radius of Ri to an outer one of R\. Along the pole, the inner and outer boundaries 
are scaled by the factor A r . Although the theory allows for fully triaxial ellipsoids, our specific applications are for axisymmetric 
ellipsoids (ranging from oblate to prolate). 
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(c.f., Messiah 1962 and the Appendix). 

As in Paper I, the scattering factors can be expressed as 



i + cos z e = - 
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and 
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so that, upon substitution in the integrals contained in Eq. (Q}, to- 
gether with the other expressions above, and using the properties 
of spherical harmonics, F sc and Z* can be rewritten as: 
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(14) 
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(16) 



Thus, the scattered flux and the Stokes parameters can be ex- 
pressed as sums of increasing order of multipole contribution. 
This helps to separate the effects of anisotropy in the flux F and 
in the density distribution n, and to study how they can lead to 
the production of polarization. 

If the functions are smooth, the summations will converge 
rapidly, so the first few terms (I and /' < 2) should provide rea- 
sonable approximations for the behavior of the polarization. Due 
to the conditions of Eqs. (0 and fllOl l, the summation over I, I', m, 
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and m' in Eqs. ( TBI and (TT3T ) will be limited. For example, if the 
stellar flux and the circumstellar density distribution of scatter- 
ing particles is symmetric about stellar and envelope polar axes 
respectively , the values of S //' mm < will be zero for odd values of 
I and 

For a spherical envelope with n(r, 9, (f>) = n(r) and an 
anisotropic light source, Eqs. ( fT4l and ( fl~5T > reduce to the forms 
discussed in Paper I. On the other hand, for an isotropic light 
source within an envelope of arbitrary shape, the expressions of 
Brown & McLean (1977) and Simmons (1982) are recovered. 

3. Ellipsoidal light source and ellipsoidal 
circumstellar envelope 

As an illustration of the preceding general formulation, we will 
take the case of a star of uniform surface brightness, as discussed 
in Paper I, but now embedded within an ellipsoidally shaped cir- 
cumstellar envelope. For an isotropic surface intensity /„, the 
flux F(r, if) can be expressed (as seen from a distant point) 
in terms of the projected area A p (i9\ <p) of the star as seen from 
direction (i?, <p). In other words we treat the star as an arbitrar- 
ily small unpolarized illuminator but describe its anisotropy in 
terms of a distorted finite shape. 

Under this approximation, the flux is given by 



F(r, &, tp) « /» AO = /, 



A p (0,¥>) 



(17) 



where /* is the isotropic intensity of the stellar surface and AO. 
is the solid angle subtended by the scattering element located 
a distance r from the star. As a concrete example, we follow 
Paper I (see their Fig. 2) by introducing star axes (a, b, c) to be 
along (X, Y, Z). From Paper I, the areal function is 



^(bcA) 2 + (acjS) 2 + (abv) 2 



(18) 



where (A,fi,v) = (cos tp sin i?, sin tp sin i?, cos i?) are the (X, Y, Z) 
direction cosines for the direction (#, tp). So we obtain 



F to (r) 



2jt 



A p (#, ip) F/Jtf, tp) sin & d-9 dip. 



(19) 



For the purpose of illustrative calculations, we now adopt 
simplifications for the description of the illumination pattern and 
envelope density profile. Following Simmons (1982), we model 
the envelope as an axisymmetric ellipsoidal shell (oblate or pro- 
late) of arbitrary thickness and constant density. The shell axis 
of rotational symmetry OZ' has inclination angle z' e with the line 
of sight (see Fig[TJ. An azimuthal angle env describes the orien- 
tation of the projected OZ' axis in the x - y plane, which is the 
plane of the sky for the observer. The density prescription for the 
envelope becomes 



ni' m '(r) = 2n(R\ -tf 2 )«o^/' Y Vm > (ie, 0e), 
where 



n(r, ■&', tp') 



no 




when r2(/j) < r < r\(p), 
otherwise. 



(20) 
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Here «o is the uniform number density of particles within the 
column bounded by r\ and ri, where 



R 



1.2 
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(22) 



and 

fi = cos if. 



(23) 



The lengths Ri and R2 are the outer and inner equatorial axis 
lengths, and A r is the ratio of the length of the equatorial width 
to the polar width (see Fig. Q]). The angle f is the angle between 
the radius vector and the axis of symmetry, which is related to 
our frames by the addition theorem of spherical harmonics. This 
explains the appearance of T/< m <(/ e , tp e ) in Eq. (l20l - see Simmons 
(1982) and Jackson (1975). Finally, is given by 



K v 
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where Pi>(x) are Legendre polynomials. 

Inserting Eqs. (fT9ll- (l24l i into Eq. ( fT6b , we obtain 

S //'mm' = I*N fl m K r T//'m'(i'e> 0eX 

where a convenient column density is introduced as 
N = 2n(Ri -R 2 )n , 



and 



./in 
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(25) 
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Note that the multipoles of the flux // m are now functions of the 
star's effective shape. The multipole coefficients of the envelope 
K/> are functions of the envelope's shape and size. Conceptually, 
the product of K/> and fi m leads to an overall net amplitude "pat- 
tern" of scattered light as determined by whether the two func- 
tions enhance or offset one another. This pattern sets the contri- 
bution level to the polarization and scattered flux as a function 
of direction about the source. The scattering phase function (in 
our case dipole scattering) along with appropriate rotations op- 
erates as a weighted filter function that acts on the amplitude 
pattern, converting it to a net observed Stokes vector signal from 
the source and envelope system. 

Describing the orientation of the star relative to the observer 
frame in terms of the Euler angles, we can choose a as zero, ft 
as the viewing inclination (/,) of the OZ-axis (i.e., the nominal 
rotation axis of the star) to the line of sight Oz, and y as the 
azimuth of the (9Z-axis from the Ox-axis as measured about the 
Oz-axis, which we denote as <p s . Thus S measures the rotational 
phase of the star relative to the observer (the same conventions 
as adopted in Paper I; see their Fig. 1). 

In the observational context, it is the normalized scattered 
flux and Stokes parameters that are usually used. These are given 
by (/so q, u) = (F sc , Q, U)/F tot , where the total flux received F tot 
comprises the combination of the scattered flux F sc plus direct 
flux from the star F*. This contribution by direct stellar light 
we denote as F* = I*A p (i s ,(p s )/D 2 . So the total observed flux 
becomes 



F ™= j\ K(4^ s ) + (0 2 //.)fsc]. 



(28) 



The general expressions for the normalized scattered flux and 
stokes parameters become 



/sc 
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Fig. 2. Plots of fio for the source flux plotted against the parameter c for the anisotropy of the source. Low c implies greater flux 
along the stellar poles, whereas large c is the opposite case. The left panel shows curves for I = and / = 2 as solid, / = 4 as dashed, 
and higher / values as dotted (unlabeled) against log c. The right panel shows a blow up around the region of c — 1 and explicitly 
labels the higher order curves. Certainly for c < 1, 1 — and 2 are most important. 
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z = 
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where t = cr(R\ - /?2)«o is the equatorial optical depth. The de- 
gree of polarization is given by p — \z*\ — \z\, and the polarization 
position angle is given by ifr = jargz. 

From Eqs. d29l ) and (|30l l, we can calculate the polariza- 
tion and the scattered flux, using the properties of the two fac- 
tors Ki> and fi m that describe the envelope geometry and source 
anisotropy. Due to the symmetry of the functions chosen to de- 
scribe the stellar flux and the scatterer density distribution, K{ 
and fi m are non-zero only for even /' and I, respectively. 

4. Applications 

As mentioned above, the multipoles of the flux /; m are non-zero 
only for even I, and the spherical harmonics for I > A are impor- 
tant mainly for fairly large angular distortions of the star from 
sphericity. As an example, Figures [2^ and b show the variation 
of fio for an oblate/prolate star distorted only along its c-axis. 
For 1 = 1' =2, the scattered flux is a sum of terms given by 
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and the Stokes parameters as 
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where 



(31) 



(32) 



(33) 



The various coefficients Ci^ , are tabulated in the Appendix in 
Tables QO 

The multipoles fi m and Kp are the determining factors for 
properties of the total scattered light. In Paper I for spherical en- 
velopes, we focused on the effect of the projected area A p (within 
the shape factor // m ), finding that the polarization increases as 
the stellar inclination z' s increases toward more edge-on viewing 
perspectives of the star. With an ellipsoidal envelope, the ratio 
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of the length of the equatorial axis to the polar axis of the ellip- 
soidal envelope A r (contained within the shape factor K/>) adds 
additional richness to the problem. For A r approaching infinity, 
the envelope will be a planar disk of scattering particles. As A r 
approaches zero, the envelope stretches to a cylindrical column. 
Note that such a form is not the same as a polar jet, since in our 
envelope parametrization the equatorial radius of the envelope is 
never zero. 

It turns out that Kq and K2 are straightforwardly derivable as 
functions of A r . There are two sets of solutions, one for prolate 
envelopes and the other for oblate ones. For prolate envelopes 
(i.e., A t < 1), the solution is 



= 2- 



3 K - 2A r 
K * = 2-TaT~ K ° 

whereas for oblate envelopes (i.e., A r > 1), the solution is 



(34) 
(35) 



K 2 = 



1 



In 



' A r + V^r 2 - 1 



3A r 

2(A r 2 - 1) 4(A? - 1) 



A r - V4 2 - 1 
1 + 2A? 



(36) 
(37) 



The Kq and K2 functions have interesting limiting behavior. For 
both prolate and oblate envelopes, the coefficients achieve their 
limiting values as A r approaches unity, with Kq — 2 and K2 = 0. 
The opposite limits are extreme distortions. As A r approaches 
zero, a prolate envelope stretches out to an infinitely long col- 
umn with Kq = 7r/2 and K2 = n/4. Note implicitly that there 
is an infinite amount of scattering mass, but it is infinitely far 
away, hence Kq and K2 remain well-behaved in this limit. As 
A r becomes large, an obate envelope becomes a thin sheet in the 
equatorial plane. The coefficients both tend toward Kq = K2 = 0. 

We next consider applications of Eqs. (fJTJ and ( 1321 to two 
particular cases. The first case deals with a scenario initially 
predicted by Owocki, Cranmer, & Gay ley (1996) for the mass- 
loss from a rapidly rotating hot star, and explored in more detail 
for the shaping of nebula surrounding Luminous Blue Variables 
(LB Vs) by Dwarkadas & Owocki (2002). For line-driving of the 
wind coupled with the effects of gravity darkening, Dwarkadas 
& Owocki (2002) predict that the equatorial wind will be less 
massive and slower than the polar flow. The second case deals 
with a a Roche lobe filling star for a binary system surrounded 
by a scattering envelope. This has relevance to short period bina- 
ries experiencing mass transfer, where one of the stars dominates 
the luminous output for the waveband of consideration. 

4.1. Rotating Stars 

For the case of a line-driven wind from a rotating star, we adopt 
the following approximations to represent the polar enhanced 
density and a reduced density from the equatorial band (c.f., 
Fig. 4 of Dwarkadas & Owocki 2002). Although the wind den- 
sity will fall off as r~ 2 or faster, our main concern is to capture 
the trend of the latitude dependence in the radial optical depth 
of scatterers, for which our ellipsoidal envelope representation 
in the form of a prolate shape is adequate. With Thomson scat- 
tering being the dominant polarigenic opacity, we introduce an 
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Fig. 3. Polarization as a function of the prolate envelope shape 
for different stellar anisotropics to represent the scenario of an 
enhanced bipolar wind flow centered on a gravity darkened star. 
The curves are: short dash for c = 1; solid for c = 0.7; and 
long dash for c = 0.4. From bottom to top, the polarizations 
are higher for more edge-on viewing inclinations (see the text). 
The bipolar wind is aligned with the c axis for the source, hence 
z' s = i e in each case shown. The optical depth in the equator is 
fixed at To,eq = 2/7T = 0.64. 



electron scattering optical depth through the equator of the en- 
velope as To,eq = "o (R\ - ^2) o~t ■ Then the optical depth along 
any other radial line is 



T 0,eq 



Vl-(l-A r V 



(38) 



where A r < 1 leads to a higher optical depth of scatterers along 
the pole than in the equator. It is worth noting that the optical 
depth averaged over solid angle is 



T = T0,e 



VT^a? 



(39) 



The ratio factor ranges from unity for A r = 1 (a spherical enve- 
lope) to n/2 for A r = (a cylindrical column). Although in this 
latter case the optical depth along the pole formally diverges, the 
average over solid angle f remains finite. 

For a rapidly rotating star, gravity darkening leads to a lower 
temperature at equatorial latitudes and polar brightening at high 
latitudes. With limb darkening included, a scattering particle at 
an arbitrary position "sees" a complicated brightness pattern at 
the star. Above the pole the radiation field as seen by the scatterer 
is centro-symmetric. In the equatorial plane, the brightness pat- 
tern of the star is banded, modulo the effect of limb darkening. 
The flux along the pole is therefore greater than in the equatorial 
plane, and so we model this trend using equations ( fTTl ) and ( fT8l ) 
with a — b and c < a as before, so the star is oblate in terms of 
its illumination. 
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We apply our theory to this scenario by adopting To jeq = 2/tt 
to obtain rough upper limits to expected polarizations in the op- 
tically thin approximation, since f < 1 for all A t < 1. Results 
are displayed in Figure [3] as a function of A, for different de- 
grees of source anisotropy as characterized by the value of c. In 
this figure the polarization q is plotted as percent polarization 
and is normalized to the optical depth in the equator To jeq - The 
short dashed lines are for an isotropic source. Solid lines are for 
c = 0.7 and long dashed ones are for c - 0.4. Each case is shown 
at four viewing inclinations of 0°, 30°, 60°, 90°. The larger polar- 
izations \q\ are for higher viewing inclinations. With the envelope 
and star axes aligned, z" s = z' e , thus the curves for the pole-on case 
have q — for all values of A r . 

With enhanced mass loss and stellar illumination along the 
poles, our model naturally predicts a polarization position an- 
gle that is orthogonal to the rotation axis of the star. Several 
LBV stars have been studied and found to have (a) net polar- 
izations and (b) axisymmetric nebulae. Schulte-Ladbeck et al. 
(1994) reported on a polarimetric study of the LBV AG Car, 
and they found strong polarimetric variability. Yet, the variable 
polarization displayed a preferred axis co-aligned with the sym- 
metry axis of the star's ring nebula, exactly opposite of what 
our model would predict for enhanced mass loss from the poles. 
A more recent study by Davies et al. (2005) arrived at a dif- 
ferent conclusion. Those authors found largely erratic changes 
in polarization, in degree and position angle, similar to what is 
observed in P Cygni (Taylor et al. 1991). In fact, Davies et al. 
found that three LBV stars had essentially random polarimetric 
variations suggesting that these winds have no preferred sym- 
metry axis but that wind clumping could explain the variability. 
On the other hand, their discussion of 77 Car is consistent with an 
enhanced bipolar wind, for which the polarization position angle 
is indeed perpendicular to the axis of the homunculus. They also 
found that R127 shows evidence for a symmetry axis. Schulte- 
Ladbeck et al. (1993) concluded the same; however, Davies et al. 
interpret the net polarization as evidence for an enhanced equa- 
torial density, not a bipolar enhancement. 

Note that although our form for r(yu) in equation d38l gives 
the correct trend with latitude about the star, the curve does not in 
general accurately reproduce the shape expected of the Owocki 
et al. (1996) model. In that paper the mass-loss rate and wind 
terminal speed scale with the latitude-dependent effective grav- 
ity, which decreases from the pole to the equator. For rotation 
speeds that are not too close to break-up, A T < 1, and the shape 
of r(/i) is a relatively good match. But for faster rotations, corre- 
sponding to lower values of A r , agreement between equation (l38l 
and an accurate treatment of the Owocki et al. bipolar winds is 
not good: the effective opening angle for the bipolar wind in our 
model tends to be too small. As a result, our treatment tends to 
overestimate the expected polarization, but it captures the overall 
qualitative trends. 

It is clear that the polarimetric behavior of LBV stars is com- 
plex. Some show evidence for axisymmetry, as one would expect 
from the Owocki et al. (1996) mechanisms, and others do not. 
They are certainly strongly variable in their polarization, likely 
a consequence of wind inhomogeniety. The lack of a symme- 
try axis for some sources may simply be the result of a range 
of rotation rates in the stars. It is also likely that the optically 
thin electron scattering may not be applicable to the high mass- 
loss rate winds of many LBVs. Large optical depths tend toward 
depolarization relative to the thin limit because of multiple scat- 
tering. The net polarization will thus be dominated by optically 
thin regions (e.g., Taylor & Cassinelli 1992). If an LBV has a 
bipolar flow that is thick to electron scattering, it may be that 
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Fig. 4. Example gM-diagrams for a prolate illuminating source 
that rotates about an axis orthogonal to its long axis. The polar- 
izations are for an oblate circumstellar envelope with symmetry 
axis coincident with the star's rotation axis. Each panel is for a 
different viewing inclination as indicated. For an edge-on view, 
u = 0, so only the ^-polarization with rotational phase is shown. 



the equatorial flow of lower optical depth could determine the 
polarization, leading to a position angle that is parallel to the 
star's rotation axis instead of perpendicular to it. More modeling 
is needed to determine whether the mechanism of Owocki et al. 
(1996) is consistent with the known polarizations of LBV winds. 

4.2. Binaries 

Binary systems provide a rich variety of scenarios involving 
distorted stars and non-spherical circumstellar and/or circumbi- 
nary envelopes. One example is /3 Lyr in which the lower mass 
companion is a late B giant star and Roche lobe overflowing 
(Harmanec 2002). In this system the more massive star is em- 
bedded in a thick disk, and there is even a jet whose origin ap- 
pears offset from either star (Harmanec et al. 1996; Hoffman, 
Nordsieck, & Fox 1998). The system also sports a significant 
kilo-Gauss level magnetic field (Skulsky 1982; Leone et al. 
2003). 

In our exploration of variable polarization in binary stars, 
we adopt fixed values of b — 2 and a — c = 1 to represent 
a highly distorted prolate-shaped Roche-lobe filling star in the 
system. Consider an envelope with an axis of symmetry that is 
parallel to the orbital axis of the binary (i.e., i e = i s ). Figure [4] 
shows variable polarizations in the case of an oblate scattering 
envelope with A r = 3 and t — 0.1. The maximum polarization 
is about 1.1%. In Paper I, this star model was considered with 
a spherically symmetric scattering envelope, and the maximum 
polarization was only 0.45%. Clearly, both the anisotropy of the 
illuminating source and the distortion of the scattering envelope 
from spherical are important for interpreting polarimetric data 
from such systems. 
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-0.4 -0.2 0.2 0.4 

q(%) 

Fig. 5. The same prolate star as in Fig. [4] but now with the oblate 
envelope not co-axial with the rotation axis. Here the envelope 
is held fixed at pole-on, and the source is allowed to tilt as indi- 
cated. 



In some cases there is evidence that the axis of reference 
for the scattering envelope is not aligned with the reference axis 
for the star. This can occur in stars with strong magnetic fields, 
such as oblique magnetic rotators (Stibbs 1950; Deutsch 1956), 
which are sometimes found in close binary systems as well. For 
a strong oblique magnetic field, the envelope's reference axis 
will be parallel to the field axis instead of the rotation axis (in 
this case for the binary orbital motion). Figures |5]-|7] shows cal- 
culations for three envelope inclinations of z" e = 0°, 45°, and 90° 
respectively. As in Figure [4] the stellar illumination is still rep- 
resented by the case a = c = 1 and b = 2 The figures show q-u 
diagrams at stellar inclinations of z' s = 0°, 30°, 60°, and 90°. The 
overall trend is for the inclination of the star to affect the degree 
of polarization more substantially than the inclination of the en- 
velope. However, the largest polarization values generally occur 
for the largest differences in the respective inclinations z' s - z e , be- 
cause the stellar flux is greatest along its smallest axes (a and c), 
while the envelope possesses a greater column density of scat- 
terers along its longest axes. 



5. Conclusions 

In this work we derive the polarization arising from an 
anisotropic point light source within an arbitrary shaped enve- 
lope. We specifically considered Thomson and Rayleigh scat- 
tering mechanisms. The mathematical analysis made use of the 
properties of spherical harmonics, which can readily be general- 
ized to more complicated cases than those discussed here (e.g., 
spots on a star). By using the first few spherical harmonics, an 
approximation was derived to represent distortions of stars and 
envelopes from spherical by varying amounts. The anisotropy 
properties of the star and envelope can be reduced to representa- 
tions as products of multipole contributions. 
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Fig. 6. As in Fig. [6] now with the envelope held fixed at an inter- 
mediate viewing inclination of i e =45°. 
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Fig. 7. As in Fig. [6] now with the envelope fixed at an edge-on 
view. 



We considered applications to rotating stars and binary sys- 
tems. For rotating stars we sought to explore the polarizations re- 
sulting from an approximate representation of enhanced bipolar 
wind flow illuminated by a rapidly rotating and gravity darkened 
star. This scenario has relevance to some LBV winds. 

We also considered a highly distorted star such as can occur 
in a binary when one of the components is Roche lobe filling 
and dominates the luminous output of the system. In this case 
the star's illuminating characteristics were modeled as a prolate 
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ellipsoid rotating about an axis perpendicular to its symmetry 
axis. Scenarios involving an oblate circumstellar envelope that 
was aligned or inclined with respect to the star's rotation axis 
were considered. 

The theory presented here is fairly general; however, there 
remain a number of improvements that should be pursued. 
Foremost is inclusion of the star's finite size, because of both oc- 
cultation and the finite depolarization effects. An initial attempt 
at this has been made by Ignace et al. (2008), who consider scat- 
tering polarization from a structured giant star chromosphere. In 
addition, our theory assumes that the radial and angular descrip- 
tions of the envelope are separable, which may not always be the 
case. Including these modifications is far from trivial but may 
be necessary in order to accurately represent complex systems. 
Certainly, the growth in the size and capability of telescopes and 
modern instrumentation, along with the increased importance of 
spectropolarimetry in ascertaining source geometries (e.g., su- 
pernova ejecta - Wang et al. 2001; Leonard et al. 2006) moti- 
vates an effort to generalize further our approach. 
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Appendix 

The Clebsh-i 

of two spherical harmonics and are given by (c.f., Messiah 1962) 



The Clebsh-Gordon coefficients Clj?£_, arise from the products 



C 



LM 



where 



a b 
a ft 



(21 + 1)(2/' + 1)(2L+ 1) 



4n 



r l 





r l 

m' M 



(-40) 



Table .3. Values of Cf mm , 
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= {-\y- b -y Va^) 

x V(a + ( fl - a V- (b + PV- (b-P)\ (c + J)\ (c - y)\ 

x \^(-\)' [t\(c-b + t + a y.(c-b + t-py 

x (a + b-c-ty.(a-t-ay.(b-t-py]}. (.41) 
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with 



A(abc) 



(a + b - c)! (b + c - a)l (c + a - b)\ 
(a + b + c+l)\ 



(.42) 



and t is an integer value for which the arguments of the factorials 
are positive or zero. The number of terms in this sum is 1 + r\, 
where r\ is the smallest of the nine numbers a ± a, b ± ft, c ± y, 
a + b - c, b + c - a, and c + a — b (c.f., Messiah 1961). 

Expression LAW is called the Racah formula and is non-zero 
under the following two conditions: 



a + ft + y = 0, 
and 

\a - b\ < c < a + b. 



(-43) 



(-44) 



Values of the Clebsch-Gordon coefficients used in this paper are 
given in Tables [71143] 
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